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Anti-charmed pentaquark from B decays 

Su Houng LeeQ Youngshin Kwon0 and Youngjoon Kwoi0 
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We explore the possibility of observing the anti-charmed pentaquark state from the ©chTr"*" decay 
of B meson produced at _B-factory experiments. We first show that the observed branching ratio of 
the to Afpir'^, as well as its open histograms, can be remarkably well explained by assuming 
that the decay proceeds first through the (or D*^) decay, whose branching ratios are known, 

and then through the subsequent decay of the virtual D° or D*° mesons to Afp, whose strength 
are calculated using previously fit hadronic parameters. We then note that the 0c can be similarly 
produced when the virtual or decay into an anti-nucleon and a 0c. Combining the present 
theoretical estimates for the ratio goNAa/O onBc ~ 13 and Qd^nBc ~ ^goNBc^ we find that the 
anti-charmed pentaquark 0c, which was predicted to be bound by several model calculations, can 
be produced via B^ —> 0ch7r’'", and be observed from the B-factory experiments through the weak 
decay of 0c —> pK^'K~■ k~ . 
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The excitement about after its first discovery by 
the LEPS Collaboration at SPring-8 0 and its subse¬ 
quent confirmation by several groups, has recently turned 
into disappointment and confusion, as increasing number 
of experiments are reporting negative results with higher 
statistics 12 0- The experimental situation is similarly 
discouraging for anti-charmed pentaquark, as the ini¬ 
tial observation by the HI collaboration at HERA|4 has 
not been confirmed by subsequent experiments 012 0 - 
While certain processes and energy ranges are ruled out, 
one can not afford to give up further refined experimen¬ 
tal search, because if a pentaquark is found, it will not 
only provide a major and unique testing ground for QCD 
dynamics at low energy, but also a basis for investigating 
many body properties of QCD at higher density. 


Given the experimental situation, one should go back 
and ask what kind of insights theoretical considerations 
can give us in the search for the 0c or the 0'''. In this re¬ 
spect, one should hrst note that the theoretical grounds 
for the heavy and light pentaquarks are quite different. 
There are on going discussions 0 0 over the validity of 
the original prediction for the mass of the 0+ based 
on the SU(3) Skyrme model [i(tI|. On the other hand, 
many theories consistently predicted that the heavy pen¬ 
taquark state is stable and lies below the DN threshold. 
The pentaquark with one heavy anti-quark was first stud¬ 
ied in Ref. |llLll2| in a quark model with color spin inter¬ 
action, with flavor spin interaction in Ei, and Skyrme 
models 00, and recently in w hich 

were motivated by the diquark-diquark|^ 0 and 
diquark-triquark|2^ picture. In the constituent quark 
model, the existence of a pentaquark state with diquark 
configuration, crucially depends on how strong the di¬ 
quark correlation is compared to the quark anti-quark 
correlation when it recombines into a meson and a nu¬ 
cleon state. Since both correlations are effectively in¬ 


versely proportional to the constituent quark masses in¬ 
volved, C Si ■ s7j HI 0 113 , the attraction is 

expected to be more effective for pentaquark state with 
a heavy anti-quark. As a simplihed example, consider 
the pentaquark picture given in ref. |2l|| . The ud diquark 
will form a color anti-triplet, isospin 0 and spin 0 state. 
Using the C determined from Ma — M^i = = 290 

MeV, one hnds an attraction of 290 MeV from the two di¬ 
quarks. This would be identical whether one has a heavy 
or light anti-quark. On the other hand, assuming the 
pentaquark recombines into a nucleon and a meson, the 
attraction expected from the diquark correlation in the 
spin 1 /2 baryon and the quark anti-quark in the spin zero 
meson would be ~ — = —430 MeV. Where C 

is determined from Mk* — Mk = — = 397 MeV. If 

rris is replaced by rric, this will become —240 MeV, from 
— Mp) = = 137 MeV. Therefore, comparing 

this to —290 MeV in the pentaquark configuration, one 
expects a bound pentaquark state only when the anti¬ 
quark is heavy. Such simple expectations are explicitly 
bourne out in the constituent quark model calculations, 
which predict a more stable pentaquark configuration 
when the antiquark becomes heavy j2d l27l 1281 l29l| . Sim¬ 
ilar results are also consistently obtained in the skryme 
model 0 and the QCD sum rule calculations . In this 
respect, the negative experimental result for the heavy 
pentaquark from the DN or D*N final state could be a 
natural consequence of its stability, and one should search 
for it from its weak decay. 

There were attempts to search for a stable heavy pen¬ 
taquark with strangeness I 23 , using high energy pion 
beam on a nuclear target. But with no realistic esti¬ 
mate on the production cross section of the pentaquark, 
it is difficult to draw any strong conclusion. Moreover, 
the QCD sum rule0 or the skyrmion approach predict 
the heavy pentaquark state to be stable only when it has 
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FIG. 1: Hadronic description for the hadronic decay —> 

AcP7r+ 


no strangenes^. 

Therefore, we propose to search for a stable heavy pen- 
taquark without strangeness, and show that the accu¬ 
mulated data at the i3-factory experiments may be suf¬ 
ficient for this purpose. To be precise, we will show, 
that the charmed pentaquark can be produced from the 
B'^ —> Oc'nTT'^ decay, and that with the most conserva¬ 
tive estimate, it can be identified through one of its weak 
decays, 0c —> pK'^'K~' k~ . 

To build up on a reliable method for estimating the 
decay, we first try to understand the baryonic decay mode 
—> A“p7r+. 

The branching ratio for this decay is well measured, 
(2.1 ±0.7) X which is the largest among the two 

and three body baryonic decays. Due to color suppres¬ 
sion factors, 7r+ is most likely produced from the W~^ 
decay, which changes the b quark to a c quark. Once 
this takes place, the remaining cu, which has the quan¬ 
tum number of a or D*°, will subsequently decay into 
A“p. In a hadronic language, this process is depicted in 
Fig.Hl Any other hadronic diagram is quite unlikely as 
the b quark has to first decay weakly for it to dissociate 
into any other baryonic modes. 

To calculate the decay rate for the process in Fig. ^ 
we first estimate the weak decay coupling GbttD from the 
branching ratio of ^ tt+D. We use the formula for 
two body decay with a form factor. 

( 1 ) 

Stt tUq 

where, |q| = ^^[{m^B “ “ 

w_d)^)] 2 . We will use similar hadronic couplings and 
form factors used previously, to fit the open charm pro¬ 
duction from p — p reaction at low energy [^. The form 
factors used are of monopole type given as, 


F{q^) 


A^ 

A2 _|_ q2 ’ 


( 2 ) 


where q^ is the D meson three-momentum in the rest 
frame of the B meson. We will take A = 300 MeV to best 
fit the open histogram to be discussed in Fig. 2. From the 
measured life time of the B^ meson of (1.671 ± 0.018) x 


10 and the branching ratio for the decay B'^ —> 

tt+D°[D*°) of (4.98 ± 0.29) x lO-^ ((4.6 ± 0.4) x lO-^) , 
we find, 


GbttD = 46 keV, 

Gb.d- = 43keV. (3) 


Note that these are just phenomenological couplings de¬ 
fined through Eq. O- 

Once this coupling is given, the three-body decay rate 
for the process B^ —> A“p7r+, represented in Fig. Q is 
given by, 




^ B+ ^AcPTr+ — 


2^Tr^m% 


dp^ 


rP D — max 

X / dpl,\M\‘^\F{pB)F{pD)\'^, 

“'Pc-min 


( 4 ) 


where 


pI 


Pd 


(m| - {pD + m^)^)(m| - {pp - m^)^) 

4to^ 

(pI) - (toa ± mpY){p\) - {m\ - mpY) 

^pI 


( 5 ) 


The range of integration is given by. 


Pd 


min,max 


{e; + eiY 

, ( 6 ) 


where 

Ep = ^{plp-ml+ml) 

^P-Kp 

EX = i^{ml-plp-ml). (7) 

‘^P'Kp 


For D intermediate state, the matrix element is given as 


|A7dP — ‘^g%pA^G'%^B 


pY - {rup ± mtYY 
(Pd - rnoY 


( 8 ) 


and for D* as, 


IMc.p = 


45i)*pAc^ 


2 

BttD* 


{2m\ + 2m^ — m 




Pd 'ml m'i 


+{m 


2 I 'm\ "mp, _ _ p'd \ 

^ ^ ^ 2 2 

- ^){{mp + tuaY -pY) 


x(r 


- TO, 


(9) 


where for simplicity, we have assumed here and in ob¬ 
taining Eq. © that the polarization sums for the D* 
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FIG. 2: Comparison of the experimental data l^ to the 
present calculation on the open histogram of the 
ArpTr"*" decay as a function of M(A7p), requiring MiAZir) > 
2.6 GeV/c^ and M{p-K+) > 1.6 GeV/c^ 



M (MeV/c") 


FIG. 3: Similar to Fig. but as a function of M{p'k'^), 
requiring M(A^7r) > 2.6 GeV/c^ and M(A^p) > 3.5 GeV/c^. 


are just proportional to 5 ^^. Using gnPAc = 13.5 and 
go-PKc = -4 and Eq. 0 , we find, 

= (2.05 + 0.51) X 10-4, (10) 

where the first(second) number come the D(D*) inter¬ 
mediate state. This lies within the experimental mea¬ 
surement of (2.1 ± 0.7) X 10-4. In Fig [3 we compare 
our calculation for the open histogram as a function of 
Pd = M{A~p) to the experimental result in ref.jl^- The 
theoretical event rates are obtained from our differential 
formula in Eq. 0 multi plie d by the total number of 152 
million B+, used in ref.|3^. The solid line is the sum of 
the contributions from D and D*. Noting that the the¬ 
oretical form factors and couplings were similar to those 
used in a totally different reaction|33j. the agreement is 
quite remarkable. In Fig. 0 we also show the histogram 
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FIG. 4: Hadronic description for the hadronic decay 

QcfLTT'^ 

in Pt^p = M(7r+p). Apart from the peak in 1.6 GeVjs^, 
which should come from the A(1600) intermediate state 
that we did not include, the other structure and magni¬ 
tude is again well reproduced. We conclude from our fit 
that the D and D* intermediate state contributions are 
important part of the baryonic decay A, and 

explains much of the detailed histogram of its decay. We 
note that such contributions have to be subtracted out 
from the data before extracting any information about 
baryon resonance mediated decay 

Let us now consider a process where a 0c can be pro¬ 
duced. The process is shown in Fig.0 We use the same 
formula as in Eq. 0 with masses replaced to those shown 
in Fig. 0 We take the me = 2800 MeV, which is slightly 
below the DN threshold, and assume it to have spin 1/2 
and positive parity. To estimate the coupling goNOa^ 
use an analogy to the 0+. 

The particle data book puts the width of 0+ to be 
0.9 ± 0.3 MeVll^. Such a small width is necessary if 
the existence of the resonance is to be consistent with 
the previous KN scattering data. Moreover, all the ex¬ 
periments reporting positive signal quote the width to 
be smaller than their experimental resolution, which are 
typically of 10 MeV. Now assuming the width of 04' to 
be 1 MeV, which is dominated by its KN decay, one 
finds that the coupling to be gKN 0 + = ijl^- Noting 
that gnpAc is estimated to be similar in magnitude to 
<7ifAfA|23i we will also take goNOc = I- Such a small 
coupling is also expected if the pentaquark wave func¬ 
tion is composed of strongly correlated diquarks with 
small spatial overlap with the DN states|2lj. More¬ 
over, we will take gDNQa/9D>NQc = goPAjgo-PA^ ~ 
1/3. With this coupling, we find the branching ratio for 
54 - _> 0 ^ 7 i 7 r 4 - to be 14.4 x IQ-^, which roughly comes 
from {goNeJgDPA.f x Tb+■ 

Once 0c is produced, and if it is unbound, it can de¬ 
cay into either D~p and D^n, and be directly observed. 
However, if it is bound, it will only be observed via 
weak decays. The dominant decay would be through 
0c —> pK'^TT~ Tr~. To estimate this branching ratio, we 
assume that it is similar to that of D~ —> Tr~, 

which has a branching ratio of (9.2 ± 0.6)%[^. 

To account for experimental acceptance and efficiency. 
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we take 70% as a rough estimate of track-finding effi¬ 
ciency including particle identification, for each charged 
particle. Then, the total efficiency to correctly find all 4 
charged tracks for a 0c decay would be (0.7)'*. Combin¬ 
ing the two i?-factory experiments. Belle and BABAR, 
we are close to accumulating 10® /B~ pairs. There¬ 
fore, the total number of expected events for 5+ ^ 
7 r+n0c and subsequently 0c ^ , would be, 

(10®)(14.4 X 10-'^)(0.092)(0.7)‘* = 32 events. (11) 

In a sense, this can be regarded as a lower limit, since 
the contribution from other possible production processes 
have been neglected and a very conservative estimate for 
goNQi: has been taken. The main uncertainty of this 
number comes from our uncertainty in the overall fit of 
Fig 2 and Fig. 3. If = 3100MeV and unbound, the 
branching ratio will only change slightly to 15.5 x lO”*", 
and its existence could be observed through strong decay 
into DN or D*N final states with event rates larger than 
that given in Eq. inj, depending on what additional final 
states will be used to identify the on shell D meson. 

Considering the fact that each R-Factory experiment 
will accumulate at least 1 x 10® R+’s every year, adding 
the already accumulated data of about 0.5 x 10®, the 
prospect of observing a charmed pentaquark state in this 
channel from analyzing the existing and upcoming data 
at each factory is quite promising. 

We have shown that the baryonic decay modes of B'^ 
can be sensibly estimated with previously determined 
hadronic parameters. Starting from such methods and 
previous estimates on the coupling of the pentaquark 
to the DN states, we find that the pentaquark can be 
produced in the baryonic decay of the R®'. Estimates 
show that in both cases, where the pentaquark is un¬ 
bound or bound, the pentaquark can be observed re¬ 
alistically through the DN or weak decay final states 
respectively, from the accumulated data at B-Factories. 
Previous searches on the charmed pentaquark were not 
combined with a realistic estimate of the cross section, 
and therefore it was not clear what to conclude even from 
a null result. Since we present a definite lower bound on 
the counting rate with which the pentaquark should be 
observed, the experimental search would be able to pro¬ 
vide a final conclusion on the existence of the charmed 
pentaquark. 
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